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We investigate the hydrodynamics of the QCD plasma using dimensionless numbers built from 
the thermodynamics and transport theory of the plasma and characteristic dimensions of the fireball 
produced in heavy-ion collisions. We find that by the usual measures, dissipation is strong, and the 
fireball is on the borderline of equilibrium. As a result, the system is richer in phenomena than ideal 
hydrodynamics would predict. One general implication is that it may be possible to get a direct 
view of the QCD plasma phase rather than having to infer its existence indirectly from signals that 
come from the freezeout isotherm after the fireball has cooled into the hadronic phase. 
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There is an emerging concordance of the first results from AdS/CFT estimates yj, lattice gauge theory [2j,|3j, and 
,— i | analysis of RHIC data |4| that the QCD plasma has small transport coefficients, i. e., small mean free times. One of 
the transport coefficients that has been measured on the lattice is the electrical conductivity Q . This is much smaller 
than weak coupling estimates would predict, implying that the resistivity is larger. Large resistivity means that if an 
external potential difference is established across the plasma, it would not drive large currents, but instead generate 
large amounts of entropy. 

In a similiar vein, if the viscosity of the fluid is small, then the implication is that Reynolds numbers are large (all 
other conditions being equal) . Large Reynolds numbers are typically associated with turbulence and other sources of 
entropy production when the fluid is subjected to external forces. On the face of it, this seems to be in contradiction 
7-H ■ to statements that the quark-gluon fluid is an "ideal liquid" , so it pays to examine this matter more closely. We do 
so, and find that the situation is complex and interesting — turbulence is unlikely to occur, but the fireball is on the 
borderline of dis-equilibrium. This conclusion is driven in equal measure by analysis of the data that the RHIC has 
produced Q, the understanding of QCD gleaned through lattice computations 0,0], and limits from AdS/CFT 0. 
More detailed conclusions are enumerated at the end of the article. 

One distinction between a liquid and a gas that seems to be easy to generalize is that in a gas the mean free path 
i-C is much longer than the interparticle spacing, whereas in a liquid the two are of the same order. The number density 
should be replaced by the entropy density, s, in a relativistic system, so one should examine the dimensionless quantity 
Q_i' called the liquidity in Q — 

43 ■ e = ™ 1/3 (i) 

where r is the mean free time. A little above the QCD cross-over temperature, T c , one may use lattice measurements 
of the two quantities on the right to estimate this dimensionless number. Using the inferred value of r = 0.3 fm from 
lattice computations at 2T C 0], and s/T 3 = 5.7 ±0.2 |[g, one finds that t — 0.8, At much higher temperatures, 

where weak coupling estimates become reliable, one has s/T 3 = 6n 2 + 0(g 2 ), and r = C(l/(? 4 logl/g), so that 
I oc 1/g 4 log(I/g) when T becomes sufficiently large. There is no phase transition between the strong coupling regime 
near T c (liquid) and weak coupling (gas) , so the nature of this "liquid" bears further investigation. 

If one takes the first order dissipative hydrodynamics of and writes it in the Eckart frame [Io|. one gets the 
Navier-Stokes' equation with a relativistic correction term. Analysis of dissipation through the Reynolds number then 
is a reasonable starting point. The Reynolds number is 

R= ^ = (e/sT)v(rT) 
j] rj/s 

where e is the energy density of the QCD matter, s its entropy density, T its temperature, v the speed of the flow, 
and r the typical size of the fireball. Since e/T 4 = 4.4 ± 0.1 for the QCD plasma at 2T C , s/T 3 = 5.7 ± 0.2 and 
v = c s « I/a/3 6], using the AdS/CFT estimate Q r)/s > 1/4tt, one finds R < 5.6(rT). In fact, on the lattice 0, 
one finds rj/s which is twice this value, giving 

R = 2.8(rT). (3) 
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This estimate is made at 2T C w 350 MeV. However, direct data on r comes from HBT analysis. At chemical freezeout 
temperature, T 175 MeV, this implies r = 7 fm. By using the freezeout values in the estimate of R, one overestimates 
it, because 77 is espected to be larger at freezeout (otherwise matter would not freeze out) and s is expected to drop, 
pushing up the value of 77/s 11]. On the other hand, if we take the temperature to be 2T Cl then r would be smaller. 
So an upper bound on R is obtained by using freezeout values of r and T in eq. ©. This is R < 18. 

If the Bjorken flow approximation is not destroyed by viscosity, then one could use the estimate T oc t _1 / Cs ~ r -3 
and r ~ r to write rt oc t~ 2 . If one neglects any special kinetics that develops at the QCD cross over, and uses 
Bjorken flow to push the estimate above back to a time when the temperature was 2T C , then one obtains R ~ 14. 
The main correction to this comes from the observation on the lattice of the slowing of sound as one approaches T c 
6]. We will later investigate the stability of Bjorken flow against dissipation. 

In non-relativistic dynamics of a liquid, one also has incompressibility. Low R and incompressibility together cause 
the motion to be reversible. This doesn't happen in heavy-ion collisions. Write 77 = ec s A where c s is the speed of 
sound and A is a typical length characteristic of momentum transport (which is c s t/3 in a non-relativistic gas). Then 



R=l±) (~)=M 



(x) <«> 



where M is the Mach number. For heavy-ion collisions, 1 < M < y/3. The lower bound comes from the fact that 
matter is expanding into a vacuum so the transverse flow essentially moves at the speed of sound. The upper bound is 
for fast particles or jets travelling through the medium. Therefore we have flow at small Reynolds number and large 
Mach number. 

Typically a large Mach number means that compression cannot be neglected and shocks are possible. A low 
Reynolds number means that energy can be dissipated during a shock. So if one has viscous shock waves in heavy-ion 
collisions, then the energy of the expansion of the fireball can be pumped back into reheating, thereby prolonging the 
life of the plasma droplet. This scenario has been studied in [L2I Il3l Il4| . 

But instead of stopping here to follow up on this, we proceed with dimensional analysis. We write the Knudscn 
number 

= (5) 
r M w 

Using the values of the Reynolds and Mach numbers presented before, we find that — A/r f«0. 032-0. 055. Another 
estimate, taking r sa 7 fm and A « 0.3 fm 0, is K « 0.04, which is consistent with this. Hydrodynamics is the limit 
of very small Knudsen number. Is the number we obtained here small in this sense? To understand this, we note that 
a finite K means that the isotherms of ideal fluid flow cannot be regarded as mathematically perfect surfaces of zero 
thickness, but are smeared out over 3-6% of the radius of the fireball. In particular, this means that the outer skin, 
from which particles can leave the fireball, is 10-20% of the volume if the fireball shape is a sphere. If the fireball 
is elongated, then the evaporative skin can be a larger fraction of the volume. Thus, in the time a low temperature 
isotherm moves half the way towards the center, almost a quarter of the material from the surface has evaporated. 
So the values of K here are not really small enough for ideal hydro, since the hot system is already pretty close to 
freezeout. 

A very enlightening way in which to consider the strength of dissipative effects is to compute the damping of sound 
waves. In the presence of viscosity, sound instensity is exponentially damped: exp(— r/T s ). The sound attenutation 
length is given by 

where the last estimate is derived from the bound conjectured in 0. At 2T C , lattice results would predict T s to be 
twice this limiting value [2J. With this estimate, one e-fold decrease in intensity would occur over a distance of less 
than 0.4 fm at 2T C . This has strong implications on the hydrodynamic history of the fireball, since both the transverse 
flow and elliptic flow are acoustic in nature, driven by pressure gradients. Such a strong dissipation of sound implies 
that the hydrodynamic prediction of the elliptic flow must be smaller. 

So, next we ask how the dissipation affects the longitudinal flow — this may be different. Using first order dissipative 
hydrodynamics 0, and doing a simple order of magnitude estimate of the relative importance of the entropy current, 
s', and the entropy s, one finds 
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Now u = 1, T = 150 MeV and r — 7 fin at freezeout. If we take the AdS/CFT lower bound, ry/s = 1/4tv, then 
s'/s s» 0.0014 near freezeout. This small number is an underestimate because rj/S has to be larger at freezeout, and 
the bound is more accurate at early times, when one has larger T and smaller r. If the entropy production were small 
at all times, one could just estimate it for Bjorken flow. Then, as before, T oc t~ 3 and r oc t, one would have s'/s oc r. 
This allows us to push the estimate above back to a time appropriate to the input value of r]/S, and shows that the 
entropy production at earlier times would have been negligible. This analysis indicates that the Bjorken solution is 
stable against first order dissipative terms even if r]/s is significantly larger than the limit conjectured in This, in 
turn, lends stability to our earlier estimate of the Reynolds number. 

While we are discussing dissipation, I make a detour to ask how big is the bulk viscosity, £. This is more or less 
unknown. In a perturbative computation it would be very small. General conditions are discussed for the vanishing 
of ( in ^5j, where a result from is quoted which says that if the trace of the stress tensor, e — 3P = f(e, n) (P is 
the pressure and n, a conserved particle number), then Q vanishes. This condition fails to hold in many cases, most 
notably when the temperature is comparable to masses of the particles or when matter contains composite particles 
and energy can be transferred from translational modes to internal excitations |l7j . In any case, a stress tensor with 
non-vanishing trace is a breakdown of conformal symmetry, and may exist in any field theory with a non-vanishing 
/3-function. In particular, in QCD near T c both these conditions may hold. 

The bulk viscosity is known for a two component model of matter where one component has very small mean free 
time compared to the other 18], and gives 



where c s is the speed of sound in the medium. Several attempts to understand the behaviour of the equation of state 
of QCD near T c have effectively involved two component models which are similiar. Moreover, at least one of the 
components in such a model is usually understood to be composite , and thereby satisfy one of the conditions for 
a large £. Since the speed of sound is now known directly from lattice computations one can form an estimate of 
this ratio. Of course, this would be correct only if QCD matter were approximately of this kind. On the other hand, if 
it is not, having some estimate available may allow us to design observables that would test or rule out this scenario. 
Using the data of in conjunction with eq. JSJ), we find that C/ 7 1 = 0.45 at T = 0.9T C . At 1.1T C , immediately above 
the transition, £/rj = 0.35. Even at 1.5T C the ratio is as large as 0.1, but it decreases rapidly above that, being totally 
negligible at 2T C and above. We end our detour on the bulk viscosity with this estimate and return to a summary of 
the evolution of the fireball. 

The viscous corrections to Bjorken flow can be appproximately absorbed into a modified stress tensor |4j|. An inter- 
esting dimensionless number to consider in this case is the ratio between the transverse and longitudinal "pressures" , 



In dissipationless Bjorken hydrodynamics, A = 1. In this form one can use the estimate P/T 4 = 1.26±0.04 along with 
the previous estimates at 2T C . Then one finds A — 1.37 if one uses the Bjorken solution to propagate the freezeout 
value of tT backwards to a time when T = 2T C . If instead one uses the freezeout value of this combination one finds 
A = 1.28. The small departure from A = 1 is another indication of the stability of the Bjorken solution. This is, of 
course, one of the bases of the analysis in Q . 

The salient qualitative results of our analysis are — 

1. The use of ideal hydrodynamics is incorrect but may not be terribly wrong in some ways. The isotherms are 
smeared out to 3-6% of the system size. However, it is bad in other ways; for example, the use of the Bjorken 
formula to extract the initial energy density is probably too optimistic. It seems possible that the initial energy 
density is significantly smaller — perhaps by 15-25%. A loss of another 15-25% by surface emission means that 
perhaps only about half of the Bjorken estimate of the energy density is available to the fireball. 

2. The transverse flow stalls. The reason is that this is a sound wave, driven by the pressure gradient between the 
core and the surface of the fireball. However, the sound attenuation length is short, since R is large. Therefore 
the pressure gradient is unable to drive a rarefaction wave efficiently. This increases the lifetime of the plasma 
over that expected from ideal hydrodynamics. 

3. This reasoning also applies to differential pressure gradients which cause elliptic flow, vi. Thus, the true elliptic 
flow must be at least 20% smaller than the ideal hydro result [2^ • In that case, the increase of V2 with s |23j 
must be re-examined to check whether it can be used to extract rj/s. 
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A = 



P + 27?/3t _ (P/sT) + (2/3)(ti/s)/tT 
P - 4t//3t ~~ (P/sT) - (4/3) (v/s)/tT ' 
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4. Bjorken flow is fairly stable under the amount of dissipation that seems to exist in the QCD plasma. The 
amount of entropy generated by purely logitudinal flow is small. However, pressure isotropy is violated by 
20-30% during Bjorken flow. Note however, that this will not drive acoustic phenomena. 

5. Any pre-equilibrium signal that one gets at the detector comes from the outer region of the fireball, and therefore 
from matter at fairly low densities. These signals should be computable in a hadronic model. 

We believe that all this is good news. If ideal hydrodynamics were a quantitatively good theory of the fireball's 
evolution, then the freezeout isotherm would have to be opaque in hadronic channels. One could not then look past 
it into the fireball. All observable hadronic quantities — single particle spectra, HBT correlations, fluctuations of 
conserved quantities, chemical abundances, would give direct information only on freezeout conditions, which lie in 
the hadronic phase. In this case, the quark-gluon plasma would be invisible; its existence would have to be deduced 
by fitting the equation of state to hydrodynamics — a notoriously ill-conditioned problem. So it is good news that the 
quark-gluon plasma is not totally invisible. One can look into the history of the plasma. However, this means that 
the tool to extract physics becomes a little complicated — one needs transport equations rather than hydrodynamics. 

I would like to thank Jean- Yves Ollitrault for discussion and a reading of the paper during a visit to CEA Saclay 
under the IFCPAR project 3104-3 called "Hot and Dense QCD Matter" . 
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